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1 Introduction 

The Q-state Potts model in two dimensions is very fertile ground for the inves- 
tigation of phase transitions and critical phenomena. For Q — 2 and 3 there is 
a second order phase transition between Q ferromagnetic ordered states and a 
disordered state. For Q = 4 the transition is also second order, but the usual 
critical behavior is modified by strong logarithmic corrections. For Q > 4 the 
transition is first order, with Q = 5 exhibiting weak first order behavior and a 
very large correlation length at the critical point. 
The Hamiltonian for the Potts model is 

W=J^(l-%, 4 )), (1) 

<i,3> 

where J is the coupling constant and < qi < Q — 1 are the Potts spin variables. 
The energies (in units of J) are evenly spaced and take on interger values in 
the range < E < Nf, where JV& is the number of bonds on the lattice. Here 
we will consider simple square lattices with periodic, cylindrical, and self-dual 
boundary conditions. 

In addition to the energy there are Q order parameters 

M, = £ (2) 

k 

which for the Ising model is simply related to the magnetization. The possible 
values of the order parameter are also intergers, < M < N s , where N s is the 
number of sites on the lattice. 

If we denote the number of states with energy E by £l(E), then the canonical 
partition function for the Q-state Potts model is 

Z Q (y)=J2n Q (E)y E , (3) 

E 

where y — e~@ J . From (3) it is clear that Z is simply a polynomial in y, and the 
analytic structure of Z is completely determined by the zeros of this polynomial, 
as first discussed by Lee and Yang[|~ 
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If we wish to study the partition function in an external held which couples 
to the order paramter, (2), then one needs to enumerate the number of states 
with fixed energy E and fixed order parameter M, Q(E,M). The partition 
function is again a polynomial given by 

Z Q (y, x) = J2 £ n Q (E, M)x M y E , (4) 

E M 

where x — e~@ h , and h is the external field. As discussed in the second^ of 
Lee and Yang's two famous papers, the zeros of the partition function for the 
Ising model in the complex- a: plane all lie on the unit circle. 

For finite systems the analyticity of Z in both x and y ensures that no zeros 
lie on the real axis. However, according to Lee and Yang, in the thermody- 
namic limit the zeros of the partition function in either the complex- a; or -y 
planes approach arbitrarily close to the real axis at the critical point, leading to 
nonanalytic behavior in the partition function. 

If the zeros lie on a one-dimensional locus in the thermodynamic limit one 
can define the density of zeros (per site) <?(0)||, || in terms of which the free 
energy per site is 

f(y) = - J g(6)log[y-y o (9)]d0. (5) 

In the critical region the singular part of the free energy is a homogeneous 
function of the reduced temperature, y — y c , from which it follows^ that g(9) 
must also be a homogeneous function for small 9 of the form 

g{9) = b { - d+v ^g{9b VT ). (6) 

This in turn implies that g{6) vanishes as 9 K as 9 goes to zero where k = 
(d — yr)/yT- On the other hand, if the system has a first order transition, k = 
and the discontinuity in the first derivative of /, the latent heat, is given by 
L = 27rg(0). Exactly parallel arguments hold in the complex- a; plane if one 
replaces the temperature exponent yx by the magnetic exponent yh- 

In a recent paper Creswick|| has shown how the numerical transfer matrix 
of Binder Q can be generalized to allow the evaluation of the density of states 
fl(E) and the restricted density of states, fl(E, M) for the Q-state Potts model. 
Similar calculations of £l(E) have been carried out by Bhanot|| in both two and 
three dimensions for the Q = 2 and Q — 3 Potts models, and Pearson Q for the 
4 3 Ising model. Beale[[lO| has used the exact solution for the partition function 
of the Ising model on finite square lattices to calculate fl (E). Bhanot's method 
is far more complex than the /J.TM and requires essentially the same computer 
resources. Pearson's method is only applicable to lattices with very few spins 
(e.g. 64) and Beale's approach makes essential use of the exact solution for the 
Ising model and so can not be used for other values of Q or in three dimensions. 
The fj,TM is quite general and the algorithm itself requires less than 100 lines 
of code. In addition, it is straightforward to generalize the fiTM to count 
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states with fixed energy and magnetization, or any other function of the Potts 
variables. 



2 Results 

The partition function for the Potts model maps into itself under the dual 
transformation 

(7) 



(8) 



where 

y- 1 

u = , 

VQ 

In the complex u-plane a subset of the zeros of the partition function tend to 
lie on a unit circle (which maps into itself under (7)); however, cylindrical and 
periodic boundary conditions are not self-dual, and this causes the zeros to move 
slightly off the unit circle. For this reason we have modified the /iTM for the 
self-dual lattice introduced by Wu et aLjnJ, so that the zeros do indeed lie on 
the unit circle and therefore are simply parameterized by the phase 0. 

In the complex x- plane the zeros of the partition function for the Ising model 
are guaranteed to lie on the unit circle by the circle theorem of Lee and YangQ , 
irrespective of the boundary conditions. 

Given that the zeros are well parameterized by a single variable, we can 
define the density of zeros for finite lattices as 

g&ftfc+l ~ AO) = - L 1 a ■ ( 9 ) 



'fcH 



In Fig.l we show the density of zeros in the complex- x plane for the Ising model 
at y — y c and y = 0.5y c - Note that at the critical temperature g tends to zero as 
one approaches the real axis, but below the critical temperature it approaches 
a constant 

2%g(0,y) = mo(y), (10) 

where mo(y) is the spontaneous magnetization. We have applied finite-size 
scaling to the density of zeros calculated in this way and find excellent agreement 
with the exact solution for the magnetization except close to the critical point 
where crossover complicates the FSS analysis. There is reason to hope that a 
more sophisticated FSS analysis will improve these results substantially. 

Finally, in Fig. 2 we show the density of zeros in the complex y-plane for the 
3-state Potts model, which is known to have a second-order transition at the 
critical point. 
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3 Conclusions 



The [iTM and its extensions offer a new way of obtaining exact information 
about finite two dimensional lattices. While the method is easily extended 
to three dimensions, memory requirements limit its use to the 2-state model 
on 4 2 x L lattices. However, Monte Carlo techniques have been developed []l2f 
which have no such limitations and show great promise in extending many of 
these results to larger lattices. Preliminary studies indicate that while most 
of the Yang-Lee zeros are very sensitive to the exact value of £l(E), the edge 
singularity and the next two nearest the critical point are not. 

In addition, it is possible to generalize the fxTM to calculate the micro- 
canonical distribution of any function of the Potts variables, and in particular 
the order parameter and correlation function. 
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Figure 1: Density of zeros in the complex rc-plane for the Ising model for L = 14, 
cylindrical boundary conditions, for y = y c (filled circles) and y = 0.5y c (open 
triangles). 
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Figure 2: Density of zeros in the complex y-planc for the 3-state Potts model 
for L = 10, self-dual boundary conditions. 
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